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Abstract 

We develop a novel and powerful method of exactly calculating various trans- 
port characteristics of waves in one-dimensional random media with (or with- 
out) coherent absorption or amplification. Using the method, we compute 
the probability densities of the reflectance and of the phase of the reflection 
coefficient, together with the localization length, of electromagnetic waves in 
sufficiently long random dielectric media. We find substantial differences be- 
tween our exact results and the previous results obtained using the random 
phase approximation (RPA). The probabilty density of the phase of the re- 
flection coefficient is highly nonuniform when either disorder or absorption 
(or amplification) is strong. The probability density of the reflectance when 
the absorption or amplification parameter is large is also quite different from 
the RPA result. We prove that the probability densities in the amplifying 
case are related to those in the absorbing case with the same magnitude of 
the imaginary part of the dielectric permeability by exact dual relationships. 
From the analysis of the average reflectance that shows a nonmonotonic de- 
pendence on the absorption or amplification parameter, we obtain a useful 
criterion for the applicability of the RPA. In the parameter regime where the 
RPA is invalid, we find the exact localization length is substantially larger 
than the RPA localization length. 
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I. INTRODUCTION 



Propagation and localization of waves of various kinds in one- dimensional disordered 
media have been studied intensively for several decades. Examples of particular interest are 
classical electromagnetic wave propagation in random dielectric media |l|fj] and quantum 
electron transport in disordered solids Much recent attention has been paid to the 

problem of wave propagation in coherently absorbing or amplifying random media 



It has been demonstrated that the wave is more strongly localized in both types of media 
than in unitary (or elastic) random media with no absorption and amplification. 

A large number of previous studies have been based on the so-called random phase 
approximation (RPA), where it is assumed that the phase of the reflected wave relative to 
that of the incident wave is uniformly and randomly distributed over all angles. One can 
easily prove that this approximation is equivalent to assuming large wave energy or weak 
disorder, and, at the same time, weak absorption or amplification. When these conditions 
are not met, an exact calculation beyond the random phase approximation is required. 

Recently, several authors have performed numerical calculation of the probability densi- 
ties of the transmittance, the reflectance, and the phase of the complex reflection coefficient 
in absorbing, amplifying and unitary random media using methods that goes beyond the 
random phase approximation |7|,|T^)-|iT[ . Their results clearly indicate that the random phase 
approximation can often lead to qualitatively wrong conclusions on the behavior of waves 
in random media. Above all, it has been demonstrated that the distribution of the phase of 
the reflection coefficient is generally nonuniform in all of absorbing, amplifying and unitary 
media. The distribution of the reflectance is also quite different from the RPA result. For 
example, the average reflectance in the RPA is a monotonically decreasing function of the 
absorption parameter, whereas the exact average reflectance is a nonmonotonic one [0. 

In the present work, we develop a novel and powerful method of exactly calculating 
various disorder-averaged quantities including the reflectance, the transmittance, and their 
probability densities in absorbing, amplifying and unitary media. We also calculate the 
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probability density of the phase of the reflection coefficient and the localization length ex- 
actly. In addition, we derive a couple of exact dual relationships that relate the probability 
densities in the absorbing case to those in the amplifying case with the same magnitude 
of the imaginary part of the dielectric permeability. From the quantitative analysis of the 
average reflectance in the absorbing case, we obtain a useful criterion for the applicability 
of the RPA. 

The outline of the paper is as follows. In the next section, we introduce the model. In 



Sec. |TJ, we describe our method of calculating the probability densities and the localization 
length in detail. We also prove the exact dual relationships that relate the probability 
densities in the amplifying case to those in the absorbing case. Results of our calculation 
are presented in Sec. [TVl Finally, we summarize the paper in Sec. |V|. 



II. MODEL 

We are interested in the propagation of a monochromatic electromagnetic wave of fre- 
quency uj in disordered dielectric media. For the sake of simplicity, we consider the one- 
dimensional case, where the dielectric permeability e varies only in one direction in space 
and the wave propagates in the same direction. We take this direction as the z-axis and 
assume the medium lies in < z < L. Then the complex amplitude of the electric field, E, 
satisfies the Helmholtz equation with the wavenumber k — u/c, where c is the speed of light 
in a vacuum, 

E(z) = 0. (1) 

The dielectric permeability e(z) is equal to 1 for z < and z > L and e(z) = l + a+5e(z)+i r y 
for < z < L, where a and 7 are assumed to be real constants and Se(z) is a real random 
function of z. We suppose Se(z) to be a Gaussian random function with zero mean and a 
white-noise spectrum: 

(5e(z)5e(z>))=g5(z-z>), (5e(z)) = 0, (2) 
4 



where (• • •) denotes an average over disorder and the constant g is a measure of the strength 
of randomness. The imaginary part of the dielectric permeability, 7, makes the medium 
absorb (7 > 0) or amplify (7 < 0) the wave without destroying its phase coherence. For 
simplicity, 7 is assumed to be uniform. The constant 1 + a is the disorder average of the 
real part of the dielectric permeability and, in general, is not equal to 1. Our method is 
applicable to the cases where a is an arbitrary nonrandom function, e.g. a periodic function, 
of z, but in this paper, we restrict our attention to the case where a is a constant. 

The above model is also relevant to the electron transport problem in disordered quasi- 
one-dimensional solids. In that case, the electron wave function ip(z) plays the same role as 
E(z) and e(z) is replaced by 1 — U(z)/E , where U(z) is the random potential experienced 
by electrons and E is the kinetic energy of incident electrons. 



III. METHOD 

We consider a plane wave of unit magnitude E(z) = ^ lk ^ L ^ z ) incident on the medium 
from the right. The quantities of main interest are the complex reflection and transmission 
coefficients r = r(L) and t = t(L) defined by the wave function outside the medium: 

e ik(L-z) + r (L) e ik(*-L) z > L, 
E{Z) = { (3) 

t(L)e~ ikz , z < 0. 



Using the so-called invariant imbedding method [FSl W\ , we derive exact differential equations 
for r and t with respect to L: 

= 2ikr(L) + % [a + z 7 + 5e(L)] [1 + r{L)f , r(L = 0) = 0, (4) 
= ikt(L) + % [a + i 1 + 6e{L)} [1 + r(L)] t{L), t{L = 0) = 1. (5) 

These stochastic differential equations can be used in calculating the disorder averages of 
various physical quantities consisted of r (= y/Re 10 ) and t (= yTe 1 ^), where the reflectance 
R = \r\ 2 and the transmittance T = \t\ 2 as well as the phases 9 and are functions of L. 
In the present work, we are mainly interested in computing the exact probability densities 



Pr(R) and Pe(0) in semi-infinite {L — > oo) absorbing and amplifying random dielectric 
media. We will also compute the exact localization length £ of the wave in absorbing random 
media defined by 

lim (logT) = -L/C. (6) 

L — >oo 

The unitary case with no absorption or amplification (7 = 0) is also of great interest. In that 
case, the probability density Pr(R) in the L — > 00 limit is clearly equal to S(R— 1) regardless 
of the strength of disorder g. But the probability density Pq{0) and the localization length 
£ depend on g in a nontrivial manner and will be studied in the present work. 



A. Probability Distribution of the Reflectance 



We obtain the probability density Pr(R) from the moments (R n ) for all integers n. An 
infinite number of coupled non-stochastic ordinary differential equations satisfied by these 
moments are obtained using Eq. (^) and Novikov's formula 2(J . It turns out that in order to 
compute (R n ) = (r n r* n ), one needs to compute the moments Z n n = (r n r* n ) for all integers 
n and h. In other words, the moments Z nfl with n = n are coupled to Z nfl with n 7^ h. 
Using the notations I = L/^o, C = k^o, a = k^a and (3 = fc£o7, where £0 = ^/gk 2 is the 
localization length with no absorption or amplification in the random phase approximation, 
we obtain the nonrandom equation satisfied by Z nfl : 

dZ n/n 



dl 



i(2C + a)(n — n) — (3(n + n) + Ann — 3n 



2 3n 2 



+ n 



— n 



~(ia-(3) -2n + 2n-l 
-(ia + /3) - 2n + 2n + 1 



Z n +l,h + n 
Z n ,h+1 ~ n 



i(iar-/3) - 2n + 2h + 1 
^(ia + ft) -2n + 2h-l 



^n— l,n 



+ nnZ n+1>h+1 + nnZ n _ x ^-i + nnZ n+x ,n-\ + nnZ n -i,n+i 

- ^n(n + l)Z n+2 ,n - -n(n - l)Z n _ 2 ,n ~ ^n(h + l)Z nfi+2 - ^n(n - l)Z n ^_ 2 , 



(7) 



with the conditions Zqo = 1 and Z n a(l = 0) = for n, n > 0. The random phase approxi- 
mation applies to the case where C ^> 1 and C ^> j3, that is 7 1. Then it is clear that Z nfl 



with n 7^ n can be neglected. This leads to the much simpler equation for Z n = Z nri = (R n ) 
obtained previously in ||: 

dZ 

( ^" -2n{n + (3)Z n + n 2 Z n+1 + n 2 Z n „^ (8) 



dl 

with the conditions Z = 1 and Z n>0 (l = 0) = 0. In this paper, we go beyond the random 
phase approximation and solve the exact equation (0) directly. 

We consider the absorbing {(3 > 0) case first. When (3 is positive, the moments Z n n with 
n, n > are coupled to one another and well-behaved for all /. Furthermore, the magnitude 
of the moment Z n n decays (more rapidly for larger (3 values) as either n or h increases. Based 
on this crucial observation, we solve the infinite number of coupled equations, Eq. (^), by a 
simple truncation method, which was developed and applied successfully to the problem of 
computing the electronic properties of quasi-one-dimensional Peierls systems by the author 
and collaborators in previous publications [pTp2f . We assume Z n n = for either n or n 



greater than some large positive integer iV and solve the finite number (= iV(iV + 2)) of 
coupled equations numerically for given values of C, a and (3. We increase the cutoff N, 
repeat a similar calculation, and then compare the newly-obtained Z nfl with the value of the 
previous step. If there is no change in the values of Z nfl within an allowed numerical error, 
we conclude that we have obtained the exact solution of Z n n. In the present work, we will 
limit our attention further and consider only the / — > oo limit, where we expect dZ nfl /dl = 0. 
Then the left-hand side of Eq. (^) vanishes and we have a set of coupled algebraic equations. 
We solve these equations by the truncation method described above to find (R n ) for every 
integer n > 0. 

It is possible to get the probability density Pr(R) from the moments by several different 
methods. We find it is efficient to use the expansion of Pr(R) in terms of the shifted Legendre 
polynomials: 

oo 

P R (R) = £ (2m + 1)(P^R))P^R), (9) 

m=0 

where P^(i?) is the shifted Legendre polynomial of order m defined over the interval < 



R < 1. The average value (P^R)) is computed using the moments (R n ) for < n < m 
and turns out to be a rapidly decreasing function of m. 

The amplifying ((3 < 0) case is a little trickier. Our method fails since \Z n a\ does not 
decay as n or h increases to large positive values. For sufficiently large values of Z, however, 
Eq. (0) is well-defined if n and n are negative. From Eq. ([j]), we can easily prove that 
Z n n{(3 = A)) = 2"-ft,- n (/3 = ~A))- Thus we obtain 

(R n )/3=f3o = (R' n )f3=-/3 , (10) 

which is equivalent to 

P R (R, (3 = -fa) = r^P R (1 (3 = A,) • (11) 

We observe that Eq. (|ll]) suggests Pr(R) to be nonzero only for 1 < R < oo in the amplifying 
case. This result does not make much sense for / 1 because in that case, we expect Pr(R) 
for (3 < as well as Pr(R) for (3 > to have a sharp peak around R = 0. Nevertheless, we 
conjecture that Eq. flTT| ) is true for all / 3> 1, that is for all L ^> £ . Once we get Pr(R) for 
(3 < 0, we can compute the moments (R n ) using the definition 

/oo 
dRR n P R (R,(3<0). (12) 

The unitary case requires a separate consideration. We expect all Z n ^s with n = n are 
equal to 1 and Pr(R, (3 = 0) = 5(R— 1) in the / — > oo limit. Therefore we have to set Z nn for 
n = N + 1 equal to 1 instead of in solving Eq. (ffi) . This changes the right-hand side of the 
last one of our iV(iV-|-2) algebraic equations, by solving which we obtain all moments of the 
form (e 4 *™-™) 6 ^. The result will be used in Sees. |III B| and [III C| in calculating the probability 



density of the phase of the reflection coefficient and the localization length exactly. 

B. Probability Distribution of the Phase of the Reflection Coefficient 

The probability density Pe{6) is most easily obtained using the Fourier series expansion 



S 



Pe(6 ) ) = 1 — V] (cos(m6>)) cos(m6>) H — V] (sin(m#)) sin(m6>). 

2?r ^=1 ' ^ m=l 



(13) 



Therefore we need the averages (cos(m#)) and (sin(m#)) for every positive integer m. 
Unfortunately, we are not aware of how to compute these averages directly when m is 
an odd integer unless (3 = 0. For even integers m = 2n, we can calculate (e 2me ) = 
(cos(2n9)) + i(sin(2n9)) using the relation (e 2md ) = Z n _ n = ((r/r*) n ). Since Z n _ n is cou- 
pled to every Z nfl with n > and n < 0, we need to calculate Y nfi = Z n -a with n, h > 0. 
The equation satisfied by Y nfl is obtained from Eq. ([?]) by a simple substitutuion: 

- 



i(2C + a)(n + n) — /3(n — n) — Ann — 3n' 2 — 3n 2 



Y 



+ n 
+ n 



-(ia- P) -2n-2n-l 
1 

~(za + /3) -2n-2n-l 



-(zck — /3) -2n-2n + l 
-(ia + /3) -2n-2n + l 



y 

■* n—l,h 
Y n ,h—\ 



- + l)F n+2 ^ - ~n(n - l)F n _ 2 ,n - ;jn(n + l)^n,n+2 - -=n{n - l)Y n ^ 2 , (14) 

with the condition F o = 1- We find that when C ^> 1, /3 or more precisely, when |2C+a| ^> 
1, P, all y n n's except Y"oo are negligible, which obviously means the phase 9 is distributed 
randomly and uniformly over the interval < 9 < 2ir. We solve Eq. fll^D in the / — > oo 
limit for (3 > using the truncation method described in the previous section. In order to 
obtain Pe{9), we also need to find (e imd ) for odd m. For this purpose, we make a conjecture 
that (e im ( e ~ n }) is a smooth function of m, the validity of which can be tested directly in the 
unitary case where we can compute (e ) for every integer m. Then the averages (e me ) for 
odd m can be obtained by a numerical interpolation between the averages for even m. Once 



we get Pe{9) for /3 > 0, it is trivial to find Pe(9) for /3 < 0. We easily see from Eq. ( ^4|) that 
^rm(/3 = A)) = ^nn(/3 = ~Po)- This implies 



Pe(0,p = /3o)=Pe{e,P = -l3 ). 



(15) 
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C. Localization Length 



In order to compute the localization length as defined by Eq. (||), we need to compute 
the average W = (logT) in the / — > oo limit. The nonrandom differential equation satisfied 
by W in the absorbing (/3 > 0) and unitary (/3 = 0) cases is obtained using Eqs. and 
Novikov's formula in a straightforward manner: 

dW 



dl 



-(1 + 0) - Re [(2 + /3 - ia)Z 10 + Z 20 ] . (16) 



In case of the random phase approximation, Z w and Z 20 vanish and the RPA localization 
length becomes || 

£rpa = (17) 

Z w and Z 2 o, however, do not vanish in general and the exact localization length is obtained 
from 

j = 1 + (3 + Re [(2 + /3 - ia)Z 10 (l -> oo) + Z 20 (l -> oo)] , (18) 



where the asymptotic values of Z 10 and Z 20 found in Sec. |111 A| are used 



IV. RESULTS 

All of our results were obtained for the cutoff N = 60. In other words, we have solved 
60 x 62 = 3720 coupled algebraic equations numerically. We have confirmed that the results 
presented in this section are exact for all practical purposes. 



A. Probability Distribution of the Reflectance 

We consider the absorbing case first. Fig. 1(a) shows the probability density of the 
reflectance Pr(R) in the large distance limit for C = 5, a = and j3 — 1, 3, 6, 10, 20. Fig. 
1(b) is the probability density of the reflectance for the same (3 values in the random phase 
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approximation (that is, C = oo), the analytical form of which was obtained previously in 



For small (3 values {(3 = 1, 3), we find that the exact Pr(R) agrees pretty well with the RPA 
probability density. As (3 increases, however, there appear remarkable differences between 
the exact and RPA probability densities. Most notably, the exact Pr(R) develops a sharp 
peak at nonzero R = R max , while the RPA probability density has a peak at R = for all 
/3 > 1. This -R max increases and the half- width of the peak decreases as j3 increases further. 
We have observed similar behavior for other values of C. In the parameter range we have 
explored in detail (1 < C, (3 < 10, a = 0), we do not find a double-peaked structure in the 
exact probability density reported in 0. 

In Fig. 2, we show the average reflectance (R) as a function of f3 for C = 1,5, 10, oo and 
a = 0. In the RPA case, (R) is a monotonically decreasing function of f3. Since 1 — (R) 
is the amount of average absorption, this means that absorption increases monotonically as 
the dimensionless absorption parameter (3 gets bigger. Though this may sound reasonable, 
it is actually a false conclusion. As pointed out in J7J], the medium with a sufficiently large 
(3 behaves as a reflector rather than as an absorber. In agreement with 0, we have found 
that (R) reaches a minimum at (3 — /3 min and is an increasing function of (3 for (3 > /? m in- 
We have checked numerically that in the (3 — > oo limit, (R) — > 1 and Pr(R) — > S(R — 1). As 
is obvious from Fig. 2, /3 min increases as C increases and can be used as a useful criterion 
for the region of validity of the random phase approximation |7],§]. That is, the RPA is 
approximately valid when \(3\ < /3 m i n (and C ^> 1 (see Sec. [IV B| )). In Fig. 3, we plot (3 m i n 
versus C for a = 0. It is a monotonically increasing function and is fitted fairly well by a 
power law function (3 min = aC b with a w 2.04 ± 0.03 and b w 0.58 ± 0.01. 

Next we show the probability density of the relectance in the amplifying case, which we 
obtain quite easily using the dual relationship Eq. (^Tj). In Fig. 4(a), we show Pr(R) for 
C = 5, a = 0, and (3 = -1, -3, -6, -10, -20. Fig. 4(b) is the RPA probability density for 



E 



P R (R, C = oo) 




for < R < 1. 



(19) 



(1-Rf 
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the same a and (3 values. In the RPA case, Pr(R) is proportional to 1/R 2 in the R — > oo 
limit for all (3 < 0, since Pr(R = 0) is a finite constant for all (3 > 0. This implies the average 
reflectance in the amplifying case is always divergent. The RPA result is wrong, however, 
because when f3 — ► — oo, the medium has to behave as a pure reflector with (R) — 1. In Fig. 
1(a), we observe that Pr(R) is finite at R = for (3 = 1, 3, 6, but goes to zero as R — > for 
sufficiently large (3 values. This observation and the dual relationship Eq. ( |TT1) ensures that 
(R) is finite for sufficiently large negative /3's. 

Finally, in Figs. 5(a) and 5(b), we show the probability density of the reflectance for 
(3 = 5, —5, a = and C — 1, 2, 4, 7, oo. In this way, we can clearly see how the exact Pr(R) 
departs from the RPA result as C decreases. 



B. Probability Distribution of the Phase of the Reflection Coefficient 

As explained in Sec. [Ill Bl we have difficulty in obtaining the exact probability distribu- 



tion of the phase of the reflection coefficient Pe(9) for arbitrary f3 values in the absorbing 
and amplifying cases. This difficulty does not arise in the unitary case (j3 = 0). Therefore 
we show the probability density in this case first in Figs. 6(a-c) for a wide range of C values 
and a = 0. For sufficiently small C"s, Pe{0) has a sharp symmetric peak located at 9 = it. 
We expect Pq(0) approaches 5(9 — n) as C — > 0. As C increases from zero, this large peak 
moves to 9 slightly bigger than tt and a small secondary peak appears at 9 smaller than tc. 
As C is increased further, the large peak keeps moving away from 9 = n and another small 
secondary peak is developed at 9 < n. The overall shape of the probability density becomes 
broader. When C ~ 0.05, the small peaks merge and are turned into a flat region. For 
C ^> 1, Pe(9) is almost constant with a small and broad peak at 9 = 57r/3 and a valley at 
9 = 7r/3. 

In the absorbing and amplifying cases, we can calculate P$(9) reliably only when 
(e* m ^~ 7r - ) ) is a smooth function of the integer m. It turns out that our interpolation method 
described in Sec. [llllj| works when \(3\ is sufficiently large compared to C or C is suffi- 
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ciently small, in other words, when the random phase approximation does not work. In 
Fig. 7, we illustrate the behavior of P$(0) for a rather small value of C (= 0.1) and a — 0, 
\/3\ = 0.1, 0.6, 1, 1.5, 2, 3. When \/3\ > C, P e {9) has a sharp peak at 9 = ir. As \/3\ decreases, 
this peak becomes lower and broader. When (3 ~ 2, the peak shifts to # < n and a new peak 
appears at 9 > ir. This new peak grows and the old peak decays as \/3\ increases further. At 
(5 = 0.1, Pe{9) is almost identical to the probability density in the unitary case (Fig. 6(b)). 

C. Localization Length 

In Fig. 8(a), we plot the inverse (dimensionless) localization length £ /£ in the unitary 
case versus log 10 C. When C >> 1, £ is close to the RPA localization length in the unitary 
case £o- 111 the opposite limit C — > 0, £ /£ g° es to zero, or equivalently £/£ diverges. Fig. 
8(b) shows the same data on a log-log plot. We note that the C <C 1 region is approximately 
linear and is fitted by a power law function £ /£ = a'C b ' with a' 0.81 and b' ps 0.65. 

Finally, in Fig. 9(a), we show the localization length in the absrbing case as a function 
of (3 for C = 5 and C = oo. We note that the exact localization length is always larger than 
the RPA localization length (= £ / (1 + P)) f° r the same (5. Fig. 9(b) shows the localization 
length as a function of C for f3 — 5. 

V. CONCLUSION 

In this paper, we have presented a novel numerical method of calculating various trans- 
port characteristics of waves in one-dimensional random media with (or without) absorption 
or amplification and used it to obtain the probability densities of the reflectance and of the 
phase of the reflection coefficient in the large distance limit, together with the localization 
length of waves. Our method is completely beyond the random phase approximation that 
has been used frequently in previous works and gives essentially exact results in the sense 
that the numerical error is unnoticeably small in all of the figures presented in this paper. 
The probability distribution of the phase of the reflection coefficient turns out to be highly 
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nonuniform when either the disorder parameter or the absorption (or amplification) param- 
eter (that is, the magnitude of the imaginary part of the dielectric permeability, |7|) is large. 
When the absorption or amplification parameter is large, the probability distribution of the 
reflectance shows behavior that is totally different from the RPA behavior. We have also 
proved a couple of exact dual relationships between the probability densities in the absorb- 
ing case and those in the amplifying case with the same |7|. The probability density of the 
phase of the reflection coefficient in the amplifying case is the same as that in the absorbing 
case with the same |7|. The probability density of the reflectance is obtained from that in 
the absorbing case with the same I7I by a simple transformation (Eq. (|TT|) ). From the quan- 
titative analysis of the average reflectance that shows a nonmonotonic dependence on the 
absorption or amplification parameter, we find a criterion for the applicability of the RPA. 
In the parameter regime where the RPA is invalid, we find the exact localization length is 
much larger than the RPA localization length. 

In the present work, we have limited our attention to the large distance limit. However, 
we could have integrated the differential equations Eqs. ([?]) and (|TJp directly to obtain the 
probability densities of the finite-size system. Our method can also be generalized in a 
straightforward manner to the calculation of the probability densities of the transmittance 
and the phase of the tranmission coefficient. Work in this direction is in progress and will 
be presented elsewhere. 
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FIGURES 




17 




FIG. 1. Probability density of the reflectance in the absorbing case for a = 0, f3 = 1, 3, 6, 10, 20 
and (a) C = 5 (exact result) or (b) C = oo (RPA). 
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FIG. 2. Average reflectance versus (3 in the absorbing case for a = and (7 = 1,5, 10, oo. 
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FIG. 3. /3 m i n versus C for a = 0. (3 m \ n is the value of (3 at which the average reflectance takes 
the minimum value. The dotted line is a numerical fit: /3 m i n = aC b with a w 2.04 ± 0.03 and 
b w 0.58 ± 0.01. 
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FIG. 4. Probability density of the reflectance in the amplifying case for a = 0, 
P = -1, -3, -6, -10, -20 and (a) C = 5 (exact result) or (b) C = oo (RPA). 
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FIG. 5. Probability density of the reflectance (a) in the absorbing case (/? = 5) and (b) in the 
amplifying case ((3 = —5) for a = and C = 1, 2, 4, 7, oo. 
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FIG. 6. Probability density of the phase of the reflection coefficient in the unitary case ((3 = 0) 
for a = and (a) C = 0.00001, 0.0001, 0.001 (b) C = 0.01, 0.1, 1 (c) C = 10, 100, 1000. 
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FIG. 7. Probability density of the phase of the reflection coefficient, P$(6), in the absorbing 
and amplifying cases for C = 0.1, a = and (a) \f3\ = 3,2,1.5 (b) \(3\ = 1,0.6,0.1. Note that 
P e (9,P)=P (9,-P). 
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FIG. 8. Inverse dimensionless localization length in the unitary case for a = (a) versus log 10 C 
and (b) versus C on a log- log plot. £o is the RPA localization length in the unitary case. 
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FIG. 9. Dimensionless localization length in the absorbing case with a = (a) versus {3 for 
C = 5, oo and (b) versus C for j3 = 5. £o is the RPA localization length in the unitary case. 
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